Behavior of colloidal particles at an air/nematic liquid crystal interface 
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We examine the behavior of spherical silica particles trapped at an air-nematic liquid crystal 
interface. When a strong normal anchoring is imposed, the beads spontaneously form various 
structures depending on their area density and the nematic thickness. Using optical tweezers, we 
determine the pair potential and explain the formation of these patterns. The energy profile is 
discussed in terms of capillary and elastic interactions. Finally, we detail the mechanisms that 
control the formation of an hexagonal lattice and analyze the role of gravity for curved interfaces. 

PACS numbers: 61.30.Jf,61.30.Hn,64.75.Xc 



Colloidal particles confined at liquid interfaces display 
rich two-dimensional (2D) phase properties [l|, |2|. The 
spontaneous formation of ordered structures such as mi- 
crocrystals has been mainly studied in simple fluids [3|-|5| 
where the self-arrangement is controlled by direct col- 
loidal interactions (electrostatic [6|, |7|, magnetic [1]...) 
and possible capillary effects. The latter might come 
from the anisotropic shape [8i] or the roughness of the 
particles [9|. It is only recently that an interest |10l - 
ll3[ has developed in the behavior of particles trapped 
at an ordered fluid interface. In bulk liquid crystals 
(LC), additional long-range interactions between parti- 
cles are present because of the partial order and elas- 
ticity. Colloidal suspensions [IJ, ll5| in a nematic ma- 
trix are thus qualitatively different from their isotropic 
analogues. They display rich self-ordering phenomena 
involving particles and topological defects. At LC inter- 
faces, complex ordered structures were also observed in 
several cases: glycerin droplets [lO|, UjJ or solid beads 
[l6| at nematic/air interfaces or microparticules at ne- 
matic/water interface [12, 13]. All these systems display 
2D hexagonal crystals that were ascribed to the compe- 
tition between a repulsion due to the bulk liquid crystal 
elasticity and a capillary attraction resulting from the 
interface distortions caused by the "nematic elastic pres- 
sure" . This new type of capillary interaction is however 
thoroughly discussed in two recent works [17|, |l8| and 
its role is not clearly established. To clarify the respec- 
tive role of the elastic and capillary force a direct force 
measurement between trapped particles coupled with a 
careful control of the LC anchoring on the beads as well 
as of the flatness of interface would be suitable. 

In this work, we present a simple technique for trap- 
ping colloids at the flat interface of an aligned thin layer 
of nematic liquid crystal. By controlling the beads den- 
sity, the interface curvature and the LC anchoring, we 
were then able to establish their respective role in the 
formation of the colloidal structures. A direct measure- 
ment of the pairwise interaction has been obtained with 



optical tweezers, which allowed us to discuss the respec- 
tive roles played by LC elasticity and capillarity. 




FIG. 1: a) Deposition of colloids at the air liquid crystal in- 
terface. The LC texture is either hybrid due to the strong 
planar anchoring on polymide and homeotropic at the air (1- 
b) or fully homeotropic on silanized glass (1-c). 

The studied systems are obtained by trapping solid 
spheres at the interface between air and a nematic LC at 
21°C. Aggregates of dry silica beads (radius R = 1.96/im 
from Bangslabs) are "exploded" by an air pulse in a 
box. The individual spheres then gently settle on a 
liquid crystal slab (flg. [1]) which avoids the presence 
of colloids in bulk. The LC layer (thickness in the 
range 10-100/im) is obtained by spin coating 4-pentyl- 
4-cyanobiphenyl (5CB from Synthon) on a glass slide 
treated with polyimide (EHC Japan) that ensures a 
strong planar anchoring. The layer exhibits a hybrid 
texture due to the strong homeotropic anchoring at the 
air interface (fig.[l]-b). Homogeneous homeotropic lay- 
ers have also been studied by using a silane treatment 
on glass with N,N-dimethyl-N-octadecyl-3amininopropyl 
trimethoxysilyl chloride (DMOAP from Aldrich) but the 
poor wettability of 5CB on silanized surfaces requires us- 
ing a surrounding glass wall as shown in Fig.[l]-c. Note 
flnally that the beads -initially dispersed in water- were 
covered with a monolayer of DMOAP following Ref. [19| , 
which ensures a strong homeotropic surface anchoring on 
5CB. They were then dried at T = 110°C before use 0. 

The colloids/LC systems were observed in transmission 
mode under a polarizing microscope (LEICA DM 2500 P) 




FIG. 2: Patterns formed by 4/im diameter colloids trapped at the air/NLC interface when a strong planar anchoring is imposed 
by the lower substrate. In thin samples, chains form along the direction of alignment. They are not observed at large thicknesses 
{d > 40/im) where the patterns strongly depend on the colloids density, going from a liquid to a loose crystal and an amorphous 
condensed state. Sketches illustrate two main possible nematic textures around the beads. 



equipped with an INSTEC hot stage (temperature reg- 
ulated at 0.1°C) and a SONY 1024x768 digital camera. 
Birefringence measurements with a Berek compensator 
were used to determine the thickness of the thinnest hy- 
brid films. The microscope also allows an accurate char- 
acterization of interfaces in reflection mode by Vertical 

thanks 



Scanning and Phase Shift Interferometries [21 



to a Mirau objective (x20) mounted on a Nano-F (MCL) 
nanopositionner focusing element. We also used optical 
tweezers based on a LEICA DMI 3000 B inverted mi- 
croscope equipped with a xlOO (NA 1.4) oil immersion 
objective, a 1064nm laser (YLM 5W from IPG Photon- 
ics) and a piezoelectric XY stage (MCL). Silica beads 
cannot be directly trapped because of the inappropriate 
index contrast in 5CB but can nevertheless be manipu- 
lated with the 'ghost' effects due to the alignment [23|. 
Tracking procedures (St Andrews Tracker [M]) were used 
to determine accurately the beads position. 

We first checked with vertical scanning interferometry 
that the beads were actually trapped at the interface. 
The top of a sphere is easily located with a Mirau objec- 
tive and the height contrast h (see Fig{T]) with the sur- 
rounding fluid additionally gives the contact angle of the 
beads at the air interface = arccos(l — /i/i?) = 31 (±2)°. 
The surrounding fluid is flat without detectable localized 
deformation (with a typical vertical resolution of a few 
nanometers). After the sample preparation, beads begin 
to organize into larger clusters, depending on their area 
density, the anchoring conditions and the LC layer thick- 
ness. The planar case is summarized in Fig.[2l Between 
crossed polarizers, a point defect close to the bead is al- 
ways observed in thin layers (d < 30/im). It is reminiscent 
of the hyperbolic defect that forms around beads of mi- 
crometer size in planar cells [1^. This defect disappears 
at large thickness, where the birefringence pattern looks 
more radial (Maltese Cross in top left). In thin layers. 



the beads spontaneously form linear chains parallel to 
the easy axis. Individual colloids are then attracted by 
those chains which grow and finally collect all surround- 
ing particles. When the thickness increases, the chains 
are much less defined and are no more observed typically 
above 40/im. At larger thicknesses, the patterns are very 
sensitive to the area density of deposited colloids. This 
evolution is shown in the top pictures. At a low colloidal 
density, a stable liquid behaviour is observed. Increas- 
ing the density typically above 1000 colloids. mm~^, crys- 
talline hexagonal domains appear and form a single crys- 
tal in a few hours. If the density is higher (above 10000 
colloids. mm~^), an increasing number of amorphous 2D 
aggregates (top right picture) are observed in coexistence 
with the crystalline structure. When the anchoring on 
the lower substrate is homeotropic the same patterns 
and density thresholds are observed indicating that the 
colloids interactions are very similar in a homeotropic 
slab and in large hybrid layers. These observations are 
somewhat reminiscent of the hexagonal lattices formed 
by glycerin droplets at the air-liquid crystal interface 
[lO|, UJj. In that case, the liquid structure at low density 
and the amorphous condensed state at large one are how- 
ever absent. Two main differences might explain these 
discrepancies. First the particles we used are solid silica 
spheres and are not deformed at the interface. Second, 
the anchoring on DMOAP is strongly homeotropic (the 
anchoring energy is W ^10~^J.m~^ uM) whereas it is 
planar degenerated on glycerin. The detailed nematic 
texture still has to be deciphered but two suggestions are 
sketched in Fig. O They are based on observations of 
beads with homeotropic anchorings dispersed inside ne- 
matic planar cells which show the presence of either a 
hyperbolic hedgehog point defect or a Saturn ring |14| . 

As said above, the nature of the interactions between 
colloids trapped at a nematic interface is still debated. 
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FIG. 3: Time dependence of the separation distance between 
two isolated particles released at 10.5 /xm and at 11/im (Inset), 
(b) Corresponding pair potential (arbitrarily fixed at ksT at 
r=20/im) derived from several beads trajectories. Inset: zoom 
around the unstable equilibrium distance Va- 



The hexagonal lattices of glycerin droplets [10| or mi- 
croparticles jl^ have been explained by the existence 
of an equilibrium inter-beads distance. The latter re- 
sults from the nematic elastic repulsion competing with a 
capillary attraction arising from the nematic pressure on 
each beads. In a recent paper, Oettel however [17] shows 
that the weak interface deformation cannot account for 
the observed effects. The possible role of many-body in- 
teractions in stabilizing the structures has also been dis- 
cussed recently in Ref. [18|. In our case, the nematic 
elasticity is clearly a key ingredient, since the hexagonal 
patterns disappear through the nematic-isotropic phase 
transition. An equilibrium distance between two parti- 
cles is however hardly compatible with the various ob- 
served structures and we focused on the bead-bead pair 
potential. With optical tweezers, two isolated beads are 
approached at an initial distance Vo and are tracked after 
trap release. We show in Fig|3l^a) a typical evolution of 
the separation distance r. As long as Tq is larger than 
Va = 10.5±0.5/im, the beads move away from each other. 
They irreversibly aggregate for Vo < Va with a final sep- 
aration of 1/im. Averaged over several trajectories, the 
interaction force fp can be obtained from the Stokes law 
25 1, and the pair potential energy Ep by its integration 
26|. The latter is given in FiglSjb). For beads with 
homeotropic anchoring, our observations therefore prove 
that an unstable equilibrium distance r^ separates a re- 
gion of attraction at short distances, and of repulsion at 
larger ones. The value of r^ is roughly constant with 
the thickness (above 40 /im) and the pair orientation, in- 
dicating that the lower planar anchoring is screened by 



the homeotropic anchoring on air. Such a pair energy 
profile differs from the one between two glycerin droplets 
expected in Ref. [10|. We first checked if the repulsive 
part between two trapped beads was compatible with a 
"pure bulk" elastic interaction in a large-distance multi- 
polar development. Whatever the exact nematic texture, 
the homeotropic anchoring at the air forbids a dipolar 
distortion [17| around a bead (for textures with cylindri- 
cal symmetry) . Dimensional analysis yields the following 
pair potential for the non-zero quadrupole moment J27j : 
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where K ^ 10~^^N is the 5CB elastic modulus in the 
one constant approximation and /3 a coefficient of order 
unity. The repulsive trajectory is then given by the com- 
petition between the drag force fy = —jv = —jr/2 and 
the driving force fp = —dEp/dr: 



r(t) = {25207TK/3^RH/-f + rl) 
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(2) 



This expression correctly fits the trajectories (Fig.[3l-a) 
with /3 = 2.1±0.2. To explain the short distance attrac- 
tion, we first examined the simple approach of Ref. [10| 
that consists in adding a capillary attraction term. The 
corresponding logarithmic dependence is however over- 
whelmed by the algebraic elastic dependence at short 
distance. A more refined approach derived from Ref. [17| 
yields the same conclusion: an additional capillary at- 
traction is unlikely responsible for the observed attrac- 
tion. A strong reorganization of the director field could 
however explain it, as already observed in the short 
distances binding observed for homeotropically-treated 
beads in bulk nematic films J28| and also supported by 
the similitude of the pair potential profile with the one 
theoretically computed for infinite parallel cylinders lo- 
cated at the nematic/isotropic interface in Ref. [29| where 
such a reorganization clearly appears at short distances. 
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FIG. 4: (a) Optical micrograph of hexagonal crystals at a 
curved nematic/air interface, (b) Geometry of the experi- 
ment (c) Measured separation distance r between colloids as 
a function of its distance x from center. 



The pairwise potential qualitatively explains the struc- 
tures formed by the trapped beads. At low densities, the 
mean distance r^ between two beads is much larger than 
Va and a liquid is observed. The hexagonal patterns spon- 
taneously form for intermediate densities where Va < Tm) 
due to the strong inter-beads repulsion. The lattice pe- 
riod is however limited by Ta and for large enough den- 
sities {rm < Ta) some particles spontaneously and irre- 
versibly aggregate. Differently from the glycerin droplets 
case, the period is therefore not univocally defined but 
changes with area density. To probe if many-body effects 
strongly infiuence this scenario, a known additional in- 
teraction was applied on the crystals. We used slightly 
convex air-NLC interfaces resulting from a non perpen- 
dicular contact angle on a cylindrical clean glass tube 
(see Fig.[T]-c). After the deposition (at a weak density), 
the beads converge to the center. After a few days they 
form stable large crystals with a period depending on 
the distance x from the center. When more beads are 
added at the periphery the system reaches a new equilib- 
rium with a decreasing mesh size until aggregates start 
to form in the center. Figure HJ-c shows the lattice pe- 
riod as a function oi x. A single amorphous aggregate 
is observed at the center. Further, the crystal appears 
with a period close to r^ that rapidly increases up to 
the crystal end. Such a variation is due to the gravity 
forces projected onto the interface that oppose the beads' 
repulsion. The resulting 2D pressure 11 = V^fp/r is de- 
rived from the inter-beads force in an hexagonal lattice 
of density n = 2/\/3r^. The balance of the forces yields: 
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where Am is the effective mass of a buoyant trapped 
bead and ax <C 1 the interface slope. Using ^ as the 
pair potential, we obtain the form 



{A - Bx^ 



-1/5 



(4) 



where A is an integration constant and B = 
Amga/7567TKl3'^R^. This expression fits the data with 
B ^ 1.1 X 10^^ m"*" (FiglU-c) while we measured a = 
92.2m~^. Although this result yields a slightly smaller 
value /3 ^ 1.4 when taking for densities pbeads ^ 2 g.cm~^ 
(provided by Bangslabs) and p^cB ^ 1 g.cm~^, we can 
conclude that the pairwise potential is enough to describe 
the observed patterns at least in a first approximation. 

In conclusion, we have exhaustively described the be- 
havior of microparticles with a homeotropic anchoring at 
a nematic/air interface and reported on the first direct 
measurements of the corresponding pairwise interaction. 
The latter is found attractive at short distances and re- 
pulsive at large ones. The long-range repulsive part is 
compatible with an elastic quadrupolar interaction. It 
satisfactory accounts for the hexagonal crystals we ob- 
served under simple or gravitational confinement. This 



same repulsive interaction could also be at the origin of 
the recently observed crystals at LC/water interface |13| . 
The origin of the attractive part is still an open question. 
The measured fiatness of the interface and the expected 
spatial dependency of the capillary interaction suggest it 
also has an elastic origin. One possible mechanism of this 
all-elastic interaction from repulsion to attraction could 
be related to a defect transformation with consequent 
colloidal binding. We are confident that this work will 
stimulate further experimental and theoretical studies to 
elucidate these fascinating interfacial phenomena. 
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